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Einstein gravitation theory can be extended by preserving its geometrical nature but changing 
the relation between curvature and energy-momentum tensors. This change accounts for radiative 
corrections, replacing the Newton gravitation constant by two running couplings which depend on 
scale and differ in the two sectors of traceless and traced tensors. The metric and curvature tensors 
in the field of the Sun, which were obtained in previous papers within a linearized approximation, 
are then calculated without this restriction. Modifications of gravitational effects on geodesies are 
then studied, allowing one to explore phenomenological consequences of extensions lying in the 
vicinity of general relativity. Some of these extended theories are able to account for the Pioneer 
anomaly while remaining compatible with tests involving the motion of planets. The PPN Ansatz 
, corresponds to peculiar extensions of general relativity which do not have the ability to meet this 

. compatibility challenge. 
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I. INTRODUCTION 



Gravity teets which have been performed in the solar system show a fairly good agreement with general relativity 
(GR) I']. The geometrical basis of GR, that is essentially the equivalence principle, is one of the most accurately 
\ verified properties of the physical world. The confrontation of observations with predictions of parametrized post- 
^ ■ Newtonian extensions of GR shows that such an extension cannot escape a small neighbourhood of GR. Tests of 
OO ' the gravity force law also result in stringent bounds on potential deviations from the standard form of GR, from 
I the millimiter range to the size of planetary orbits . There are however a few indications which may be seen as 
challenging GR. 

Observations at galactic and cosmological scales have led to introduce large amounts of dark matter first, of dark 
energy then, in order to fit the predictions of GR. As long as these dark components are not seen by other means, the 
' observations can as well be attributed to modifications of the gravitational force law at very large scales @,|l,|a,H0il3l- 
This idea can be tested by following the motions of deep space probes reachingthe outer part of the solar system 
and the first probes of this kind have in fact shown anomalous accelerations 0, |0 . The anomaly has been observed 
f : as unexplained residuals of the Doppler tracking data from Pioneer 10/11 probes with respect to predictions of GR. 
They are usually described as an anomalous acceleration directed towards the Sun with an approximately constant 
amplitude ~ 0.8 nm s^^ over a large range, 20 to 70 astronomical units (AU) , of heliocentric distances. Up to now, the 
^ anomaly has escaped all attempts of explanation as a systematic effect generated by the spacecraft or its environment 
k>( I [lllIT^ . This status has motivated a lot of studies aimed at a better understanding of deep space navigation as well as 
''^ • investigating the possibility that the Pioneer anomaly be the first hint of modified gravity laws at large scales 0, j 
■ and led to proposals for new missions [Tsf. 

A new kind of extensions of GR has recently been proposed with the remarkable property that they are able to 
account for the Pioneer anomaly while remaining compatible with other gravity tests performed in the solar system 
[T^IT^ . This feature has been analyzed in the context of a linearized treatment of gravitation fields, which restricts 
its appHcabihty. It is the aim of the present paper to overcome this restriction by studying the new framework without 
assuming such a linearization. Before going into the heart of the subject, it is worth recalling the motivations for 
a generalization of general relativity, which are rooted in its now long standing confrontation with Quantum Field 
Theory (see more references in p^flTj '). 

When considered as a field theory, GR appears to be a good effective theory at intermediate scales where accurate 
observations or experiments are available. However, radiative corrections due to the coupling of g ravity field to the 
other fields naturally lead to consider GR as immersed in a family of fourth order theories [Tslligl lioj with different 
behaviours at shorter or longer scales. This larger family shows the advantage of exhibiting renormalizability as 
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well as asymptotical freedom at high energies |21|, 123. |23| . This suggests that extensions of GR can in principle 
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be tackled by using the powerful tools of modern renormalization group equation • A- complete solution 

of renormalization group trajectories appears as a formidable task at the moment but preliminary studies already 
confirm that GR can be a considered as the limit at intermediate ranges of a well defined field theory with different 
behaviours at shorter and longer ranges. As this research program is still far from completion, an alternative approach 
is to explore potential extensions of GR from a phenomenological point of view. 

This view meets the aims of the parametrized post-Newtonian (PPN) formalism 0: theories lying in the vicinity 
of GR may indeed be explored through the anomalies they can produce in gravitational observations or experiments. 
Of course, these theories have to obey an important criterium of compatibility with the stringent bounds set by the 
variety of tests already performed in the solar system. As we shall see in the following, PPN models are peculiar cases 
of the more general extensions studied in the present paper. They cannot account for the Pioneer anomaly while 
remaining compatible with planetary tests ^3 1 in contrast to the general extensions of GR studied here [T^ [l^ . 

These extensions can be characterized as 'post-Einsteinian' for the two following reasons. First, the core of Ein- 
steinian theory of gravity is fully preserved : the gravitational field is identified with a metric tensor g^i, and motions 
are described by associated geodesies. As a consequence of this geometrical foundation of the theory, the weak equiv- 
alence principle is unaffected. As is well known, this principle is well verified at distances ranging from the millimeter 
in laboratory experiments to the sizes of Earth- Moon __28J or Sun-Mars orbits |2^|33. The relative accuracy, 
better than 10^^^ for some of these tests, is g ood enough to discard any interpretation of the Pioneer anomaly from a 
violation of the equivalence principle pXl l3ll | . This does not mean that the equivalence principle is an exact property 
of nature, but only that its potential violations js^ls^ are smaller than the deviations from GR studied in the present 
paper. 

Then, the coupling between curvature and energy- momentum tensors differs from its standard form in GR [34j . 
More precisely, the Einstein tensor £"^1^ = R^^ — g^^R and the stress tensor T^i, are simply proportional to each 
other in GR whereas they are related through a more general linear response in the new 'post-Einsteinian' theory. In 
other words, the Newton gravitation constant Gn which expresses the proportionality of E^i, to T^^ in GR is now 
replaced by two running coupling constants, which depend on scale and differ in the sectors of traceless and traced 
tensors. Note that and T^^ are transverse tensors with their linear coupling preserving this transversality '35|. 
The consequences of this modification have been studied in |l6lll7j | within a linearized approximation, but they will be 
investigated in the present paper without such a restriction. The gravitational effects on orbits and light propagation 
will be calculated with the non-linearity of the theory fully taken into account. 

As a result of such a study, we will show that the 'post-Einsteinian' extensions of GR show the ability to account 
for the Pioneer anomaly while remaining compatible with the gravity tests performed on planetary orbits. The 
phenomenological discussion will be presented in terms of two potentials representing the two running coupling 
constants in the space-time domain as in [T^ ll^ . The first potential, which generalizes the usual Newton gravitation 
potential, will have to remain very close to its standard form in order to preserve the good agreement between GR and 
gravity tests performed on planetary orbits, either by checking out the third Kepler law, or by looking at the precession 
of the perihelion of orbits jSg. Meanwhile, the second potential will open additional freedom and offer the possibility 
to accomodate a Pioneer-like anomaly for probes with eccentric trajectories. It will also lead to other observable 
anomalies, in particular in time delay, deflection or Doppler tracking experiments on electromagnetic sources passing 
behind the Sun |33,|3i|3ili3■ 
In the next sections we will present the modified gravitation equations (sect. 2) and their solutions in terms of metric 
fields (sect. 3), using the fact that the deviation from GR must remain small in order to preserve compatibilily with 
solar system tests of gravity. We will express the small variation of the post-Einsteinian theory from GR in different 
manners, in particular two potentials representing the modifications of the metric fields, two independent components 
of Einstein curvature tensor which may now differ from zero outside the gravity source, or two running gravitation 
constants exhibiting a scale dependence. The relations between these equivalent representations will be written down 
explicitly. 

We will then study some potentially observable consequences of this small variation of the post-Einsteinian theory 
from GR, looking at geodesies followed by probe masses (sect. 4) as well as electromagnetic ranging between Earth 
and such probes (sect. 5). The Pioneer-like anomaly, which has been predicted to exist for probes having eccentric 
motions in the outer solar system [T^ [l7j , will be confirmed here without the restrictions of a linearized approach. 
Perturbations of bound geodesies will be evaluated (sect. 6), and the absence of observation of anomalies of perihelion 
precession with respect to GR will be translated in terms of constraints on the deviations from GR met in the solar 
system. Observational constraints will finally be gathered (sect. 7), leading to the conclusion that the theoretical 
arguments presented in the beginning of the paper and the observational ones discussed then, when taken together, 
appear to favor a simple extension of GR introduced in |l6l Il7| as accounting for the Pioneer- like anomaly. 
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II. MODIFIED GRAVITATION EQUATIONS 



As discussed in the introduction, we use the high accuracy of tests of the principle of equivalence to consider only 
metric extensions of GR. Moreover, we focus our attention on a static and isotropic metric representing space-time 
around a punctual stationary source. In the solar system, this assumption amounts to disregard the effects of rotation 
and non sphericity of the Sun which affect significantly inner motions [4ll |43 | but have a small influence on outer 
ones. This assumption notably simplifies the description since metric fields only depend on two functions of one 
variable, the radius coordinate. This dependence will be written below in different coordinate systems, in particular 
Schwartzschild or Eddington isotro pic coordinates commonly used to represent the Schwartzschild solution Q| or 
its PPN extensions jil El (ii • 

In this section, we write down these equivalent representations of the metric fields, together with their relations to 
curvatures. We then deduce the form of the equations replacing those of GR, introducing anomalous gravitational con- 
stants in the two sectors of traceless and traced tensors which have the same phenomenological content as anomalous 
parts of the metric. Note that these relations, obtained in the static isotropic case, constitute a first approximation of 
actual anomalies. They will have to be checked out in future works by performing more complete calculations taking 
into account perturbations of the static isotropic case, due to the structure and rotation of the Sun as well as to the 
presence of planets (see as an example the detailed calculation of the Pioneer anomaly within the context of GR in 

. . . . 

A static isotropic metric may be characterized by two functions, corresponding to time and radial components 
written in terms of isotropic coordinates 

ds2 ^ g^Q^r^dt^ + grr{r) {dr^ + r^idO^ + sin^edip^)) 

dng^iu = de9t_,v = <9^.9piy = (1) 

(co ordinates x = {ct, r cos cos ip, r cos 9 sin ip, r sin 9) with c the velocity of light). The conventions are the same as in 
|43| with some definitions and general properties gathered in^ An equivalent representation of the same spacetime 
may be given in terms of Schwartzschild coordinates 51] 

ds2 ^ g^^[f)c^dt^ + grr{r)df^ - {d9^ + sm'^9dip'^) 

dng^iu = de9t_,v = ^^^.g^,^ = (2) 

The metric tensors, g^^ for the first coordinate system and for the second one, do not depend on the variables 
t, 9, ip which are the same in both systems. The relation between radial coordinates r and r depends on the metric 
itself 

dr i—^df ( rdrgrr\^ df 

— = V-9rr— = 1 + — (3) 

r r \ Zgrr / r 

Equations ||2Jl also provide the functional relation between the two different representations {goo^grr) and (goo^grr) of 
the degrees of freedom of gravitational field. 

General expressions of the curvature tensors are given in^as well as Bianchi identities between these expressions. As 
the static isotropic metric from which they are calculated, these curvatures contain only two independent components. 
These two components may be chosen for example as the scalar curvature R and the single independent component 
W contained in the Weyl tensor 

W = W'^^r = W^'^e , Woe = W°^o^ = = W^'^r^ = - y (4) 

These two quantities are derived from the metric components, for example in terms of isotropic coordinates, 

i?(r) = -—[fdr + -]dr (21n|(7,,| + Ingoo) 

9rr \\ rj 



1 / 9 2 

+ - {^{drlnlgrriy + drln\grr\dr\ngoo + (drlngoo" 



W(r) — ( dr ~ drW ] drW, w = -In 



500 



(5) 



The advantage of this representation is that the scalar curvature R and Weyl amplitude W behave as scalars under 
reparametrizations of the radial coordinate 

R{r) = R{f) , W{r) = W{f) (6) 
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One may use equivalently another representation, with the two components Eq and E"^ of Einstein tensor as indepen- 
dent degrees of freedom for gravitational field, 

E^{r) = —{fdr + -) a,.ln|g,,,| + \ (9.1n|.9,,|)' 

9rr \\ r J 4 



9' 



rr 



Kir) = — ( - idr\n\grr \ + drlngao) + \dr\n\grr \ (9rln|5„| + 25rln5oo) 

grr \r 4 

^ V 9rr .900 / 

Note that the remaining components of the Einstein tensor Eg and E^, or Eg and i?^, are equal and determined by 
Bianchi identity 

V,,E^ = drEi: + i (a.lngoo) (E^ - E°) + (dMr^grA) {e; - - (8) 

We now write generalized gravitation equations as a linear response relation between the Einstein curvature tensor 
and the energy- momentum tensor |16i . il7i | 

E^^{x) ^ I d^x' {x)t{i^,x')T^{x') (9) 

General relativity corresponds to the simple Einstein-Hilbert relation 34], that is also to a local expression for the 
linear response function (Gjv is the Newton gravitation constant), 

^ '-^Tt^ (10) 

[ix)7xi^,^X, ^^^{5i5{ + 5'i6i)5^^\x-x') 

Now the linear response equation, arising for example from radiative corrections of GR through the coupling of gravity 
to quantum fields js^ . should contain a non local part, equivalent to a scale dependence of the response, 

E':; ^ + 6E'i , SE^Ax) ^ j iSx)Zi^, ^')T^i^') (H) 

The integral equation then takes a simple form when the energy-momentum tensor of the source is static and punctual 
(M is the mass localized at the origin) 

T>^{x)=5'XMcH^^\^) 

SEi^ix) ^ Mc' f cdt' {5x)tlix. x') , x' ^ (cf, 0) (12) 



In preceding equations, [E]^, denotes the standard GR solution and 5E the deviation of the post-Einsteinian solution 
from the standard one. At this point, we must emphasize that observations in the solar system are compatible with 
GR, with the notable exception of Pioneer anomaly, and that deviations from GR are bound to remain very small 
0. We may also note that theoretical reasons point to the same conclusion, namely to restrict the attention to 
extensions remaining close to GR. As a matter of fact, radiative corrections show a better behaviour with respect to 
renormalization when computed in the vicinity of a solution of GR rather than in the vicinity of Minkowski metric 
|2lj| . Fourth order theories surrounding GR also show asymptotic freedom at high frequencies l23ll and preliminary 
studies of renormalization group trajectories, which are restricted to a single sector up to now |24l. l2,'l l2fi|| . suggest 
that the full theory remains close to GR over a large range of length scales, precisely those intermediate scales which 
are tested in the solar system . 

For these reasons, we will only consider extensions lying in the vicinity of GR. This entails that the anomalous terms 
5E or (5x)i well as the anomalies to be observed in gravitational phenomena, are small and may be calculated as 
perturbations around GR. As the latter corresponds to Einstein curvatures vanishing everywhere except on gravity 
sources, that is the Sun in the solar system, this assumption allows one to simplify the Bianchi identities (see the [X|l . 
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Further simplifications arise when using the fact that potentials or curvatures have relatively small amplitudes in the 
solar system (see the^J- For example, the transversality of Einstein tensor then takes its simple Minkowskian form 
and the tensor can be decomposed on the sectors of traceless and traced tensors as in the linearized approximation 
|l6l IT^ . The generalized gravitation equations may then be written in terms of two gravitation constants which differ 
in the two sectors and depend on momentum or length scale. Two running coupling constants G° = Gn + 5G'^ and 
— Gn +5G^ thus replace the unique Newton gravitation constant Gtv, introducing a dependence on the conformal 
weight, different in the two sectors, as well as on the length scale. The information contained in the anomalous 
running coupling constants 6G^ and 6G^ is equivalently encoded in the spatial variation of the anomalous Weyl and 
scalar curvatures 5W and 5R of a solution or, also, of its two anomalous components SE'^ and SE"!^. 

These equivalent formulations of the extensions of GR will be used in the following. General relations between them, 
assuming only the idea of perturbation in the vicinity of GR, are to be found in the next section. These expressions, 
obtained around GR, preserve their full non linear dependence in Newton gravitation constant Gm- Simpler relations 
are obtained in the linearized regime, which corresponds to a weak Newton potential or, equivalently, to a perturbation 
around Minkowski spacetime. An explicit correspondence may be deduced in this case between the running coupling 
constants 5G^ and 5G^ and anomalies in the solutions (see the El- In particular, these anomalies can be written in 
terms of two potentials which were introduced in the linearized approach of |16l Il7l | but can be used also in the non 
linear approach of the present paper (see the next section) . 



III. POST-EINSTEINIAN METRIC 



In this section, we discuss the post-Einsteinian metrics which represent the spacetime solution for extensions of 
gravitation theory in the vicinity of GR. These extensions are best characterized by the anomalous Einstein curvatures 
5E'^ which may now differ from zero apart from the gravity source. We express metric tensors in both systems of 
Schwartzschild or isotropic coordinates. As emphasized in the Introduction, we fully account for the non linearity of 
the relations between curvature and metric tensors. 

We begin with the solution of the non perturbed theory, that is to say GR, 

[£;;]^^ = 8..5%^5(3)(x) , (13) 

The constant k is related to the Schwartzschild radius and the metric can be written in terms of Schwartzschild or 
isotropic coordinates (cqsl^lorQ]) with 

[ffoolst = 1 - 2k?2 = -— ^ , u=\ (14) 

[9rr\st 

1 - KU/2V _ 1 



Metric components have been defined in terms of a single Newton potential ku or nu, however expressed in terms 
of different radial coordinates r or f in the two coordinate systems; inverse radial coordinates u and u have been 
introduced. All curvatures but Weyl one vanish outside the source for this GR solution 

(1 + K7i/2) 

Extensions of GR are well characterized by the fact that Einstein curvatures no longer vanish apart from the gravity 
source. The anomalous non zero values 5E'^ are entirely characterized by two independent components, for example 
5Eq and SE^., the other components being given by the Bianchi identity ((HJ which is now read 

drSEl + i [a.lngool^t {^1 ~ SE°) + [dMr^9rr\] {SE^ ~ 6EI) = (16) 

5Eq and (5i?^ can be seen as representing the anomalous running coupling constants (sec the previous section). They 
may be written in both coordinate systems, taking into account the fact that inverse radial coordinates (u, u) satisfy 
relation 1)14(1 in GR, but a modified relation in extended theory 

^ Sgrr , , 

— m = 777 — r- (17) 
u 2 [5rrJ^^ 
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Anomalous metric components and Einstein curvatures are equivalent representations of the extended theories. 
Their relation may be rewritten, using Schwartzschild coordinates, 

goo = [.90o]st + '^500 , 9rr — [.9rr]st + ^9rr 

-r^du'+ / ,3'^ , 18 
u'^ J u-^ [5oo]st 

For the forthcoming discussions, it will be convenient to introduce two potentials 5^i^ and (5$p through the definitions 

2u 
(1 ^2ku) 



5grr _ 


u 




[500 1st . 


Sgoo 


f du 




[300] ft 



^500 - 2 1 - 2ku) / „ du 19 

J (1 - 2k.uY 

Let us also note the following relation, which will appear useful in the following, 

5g'^Q = 2(S$^ + 2K5(5oo5rr) (20) 



In the linearized approximation, all corrections in ku may be disregarded and the simple relations used in |16l Il7| 
are recovered. Clearly the full equations 119|) describe non linear effects of the Newton potential and have therefore 
a larger domain of validity than their linearized expression. The precise form of the non linear version (|19|) has been 
chosen so that they lead to simple expressions for the anomalous Einstein curvatures 

SE; = 2u^S^'p (21) 
Finally similar relations may be written in terms of isotropic coordinates 

-2(1-2ku)2 / ^ _ ^ du 



[ffrrlst J {l~2Ku)i 

Sgoo ^ f S^^-2num,,^ 2^ n^i.-^y,^ (^2) 



[.goolst J (1-2ku)2 {1-2ku)2 J (1-2ku) 



or, alternatively. 



Sgrr r, -\i /" du SEq 



[grrU (1 - 2k.u)2 J m'" 



du' 



loo 



[5. 



'00 



St 




(1 - 2ku)2 J {I- 2ku) 



(23) 



To sum up the results of this section, relations (|19|) and \i'6\ describe the looked for correspondance between 
different representations of extensions of GR: anomalous Einstein curvatures and Schwartzschild metric are given in 
terms of potentials 5^n and (5$p by relations which take into account the non linear effects of standard Newton 
potential. This feature will be preserved in the following sections, when consequences will be derived for geodesic 
motions and light propagation. At this st age, it is worth comparing generalizations obtained in this manner with the 
widely used PPN Ansatz [H IH IHIil [H llOl • This PPN Ansatz is written inO as a particular case of the more 
general extension presented here. Non zero values are found for Einstein curvature 5E^ apart from the source, at 
first order in k, but not for 5Eq which vanishes at the same order. It follows that the PPN Ansatz does not have the 
ability to account for the Pioneer anomaly, since the latter is associated to a curvature 5Eq which does not vanish. 
Moreover, the scaling laws are different for the PPN Ansatz and for the extension accounting for the Pioneer anomaly 
(see sect. 7). 
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IV. MODIFIED GEODESICS 



We derive in this section the geometric elements associated with post-Einsteinian metric which wiU be needed for 
discussing gravity tests. As in the two previous sections, we consider extensions lying in the vicinity of GR so that 
any expression will be obtained as the sum of a standard GR value and of a small anomaly, to be treated as a first 
order perturbation around GR. 

As in any metric theory, motions obey the geodesic equation 

The covariant derivative D is determined by the Christoffel symbols associated with the metric or (|21l (see 
the El for definitions). Accelerations undergone along the geodesies may be defined as Christoffel symbols projected 
along the motion 

9fJ.^-^ = -^tJ^^i^pu'^u'' = -r^^ss (25) 

Each geodesic is contained in a constant plane which may be chosen as = ^. A given geodesic is then determined 
by two constants of motion, the energy E and angular momentum J, for a test mass m 

dt ^dtp , 

mc-'goo— = E , mcgrrr -j- = J (26) 
ds ds 

Accelerations are then written, using isotropic coordinates, 

2 



1 / / , f dr\ „ / 900 



drgoO - dr{g009rr) 3- -Or 



2goo \ V"*/ \T yrr J m'^c'^ I 

r,,. = p^^W (27) 

\r grr / mc ds 

They may also be expressed as sums of a standard GR expression and an anomaly, written at first order in the 
deviation from GR, 

XT f ^900 

2 V V 5ooLt/ 5ooLt / Vrf« 



'^.900 _ _S9rr_\ _ / [goo] St \ g\ J"^ 



2 



^r... = dr ( (28) 



[grr\stJ rncds 

Note that inserting the explicit form of metric perturbations H23|l leads to expressions of perturbations of geodesies 
in terms of the anomalous Einstein curvatures. 

In order to discuss electromagnetic ranging, we will also need to express the perturbations of time delays. The 
latter are represented by a two-point function T describing the elapsed time along the propapagation of a lightlike 
signal from a spatial point (ri,(^i) to another one (r2,(p2). Equivalently, this function represents the propagation in 
a static isotropic metric of a null geodesies taking place in the plane ^ = f . This problem has been solved in its full 
generality in |17| and we reproduce here the solution in terms of integrals along null geodesies 

Pdr/r'^ 2w{r)^_9rr. 



/g2iu(r) _ £1 500 



1-2 „2w{r)^^ 



'7'i2 (j-i, Vi;^'2, V2) = / — . (29) 
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The two endpoints of the null geodesic have been supposed to lie on the same side of the gravity source. When this is 
not the case, the time delay can be written as the sum of the delays between each endpoint and the point of closest 
approach. In all cases, the latter is determined in a self-consistent manner from the two endpoints by H29(l . Note 
that the integrals H29|) only depend on the ratio (—grr/gon) as a consequence of conformal invariance of null geodesies 
propagation. 

Expression H29(l of the time delay function may then be expanded as the sum of the standard GR expression and 
of an anomaly calculated at first order 



cTi2 = c [Tl2]^^ + cSTi2 , cSTu = / 5w{r)cd [T]^^ 



e 



L 1st 



cd[TU^ _^ , ^^=o r-^-TI^ (30) 



c [T]^^ is the standard optical distance in GR and cST the first order variation in the extended theory. This deviation 
only depends on the perturbation 6w of the conformally invariant ratio (—grr/goo)- The relation l|3U|) translates into 
quantitative terms the qualitative idea of a dephasing of phasefronts by the anomalous Weyl curvature. 



V. PIONEER-LIKE ANOMALY 



We come now to an essential property of the new framework, namely the prediction of a Pioneer-like anomaly for 
probes following eccentric trajectories in the outer solar system. For that purpose, we first recall the results obtained 
in ^3 for writing the effects of the modified metric on Doppler shifts. We do not repeat the part of calculations which 
already included non linear effects of Newton potential. 

Pioneer probes were followed through a Doppler velocity built up on radio signals exchanged with a station on Earth: 
an up-link radio signal emitted by the station at Xi — (rf^,xj~) and received by the probe at X2 = (ct2,X2) was 
instantaneously transponded to a down- link radio signal then received by the station at = (ct]*", x]*"). Both Earth 
and probe followed geodesies, with respective velocities uf,U2. The gauge invariant Doppler observable, denoted y, 
is defined as the ratio of clock rates of a reference clock at the departure and arrival of the radio signal [13 

= ^ , ds'' = g^.dxt'dx'' , u'^ = ^ (31) 
dsi ds 

The Pioneer effect is observed when looking at the Doppler shift H31|l as a function of time. As in we write it as 
an acceleration 



2ads = -c^dy (32) 

The definition of a depends on the choice of the reference clock delivering ds whereas the definition of ads does not 
depend on this choice. 

All these observables may be determined from the time delay T12 on up- and down-links and on the motions of the 
endpoints. Motions of Earth and probe have vanishing covariant accelerations Dui — Du^ — Du2 — 0, corrections 
associated with the position of stations on Earth being taken into account separately (these small corrections are 
supposed to be linearly superposed to the anomaly and to disappear in the end of the calculations). With these 
assumptions, the following exact relation is obtained [T^ 

2ads _ (wj" ■ Dn^) ^ [u^ ■ Dr4) _ (^2 ■ Dn2) _^ {u2 ■ J'n^) ^^^^ 
{u^ -n^) (u^-n^) (u2 • nj) {u2 ■ n^) 

represent the directions of wavevectors at endpoints and • denotes a 4-dimensional scalar product 

{nf)o = -(nJ)o 
(u ■ n) = u^n^ 

The Doppler acceleration (|33|l has a twofold dependence versus the metric. A first dependence is due to the 
Christoffel symbols which enter the covariant derivatives and a second one to the variation of the time delays which 



dT^2 



(34) 
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enter expressions (|34() . When evaluating the Pioneer anomaly, that is to say the variation of a due to the metric 
perturbation with respect to GR, we thus obtain the anomaly as a sum of two terms 

6a = Sax + Sar (35) 

Sa-T corresponds to the variation of the time delays on radio links whereas Sar enters through the variation of the 
probe trajectory 



26aTds2 ^1 1 cdST{x^ ,X2) ^ 1 cdS7'{x^,X2 
(it2 • n^) ds2 (u2 • n^) ds2 



2SaTds2 _ , - (ST^^ • n J ^ (5r+ 



n 



("i ■ rii ) « • ) 



(U2 • ) (u2 • ri2 ) 

We have specified the reference clock to be on the probe ds = ds2 , dsf being associated with the clock on the station 
orbit. Calculating at first order in the deviation from GR, we have replaced all quantities but 6T and 6T by their 
zeroth order approximation. 

From now on, we furthermore specify that the Pioneer probes lie and move in the radial direction. Using (|27() . one 
sees that the contributions of the Christoffel symbols for the Earth may be neglected. The anomalous accelerations 
(|35|l may thus be rewritten 

— IT- — -7-n (cc?(5T(x7,X2) + cd(5T(xt,X2)) 



26aT 



Inserting explicit expressions of time delays (|30|l and Christoffel symbols (|28|) . we rewrite these two contributions in 
terms of metric perturbations (functions are evaluated at the probe radial coordinate r2) 



2San 



2SaT 




Or 7 T- [grrlt + T^^r + OrOgrr -j— 

V 5ooLt/ 5ooLt / \ds2 



(38) 



Using expressions H22|l . we may finally obtain the anomalous accelerations in terms of the two potentials 5^m and 
(5$p. This leads to the non linear generalization of the expression obtained for the Pioneer anomaly in the linearized 
approximation |l7j| . This expression is recovered in the particular case of probes located in the outer solar system, 
where Newton potential is small. We will use it in the end of this paper for discussing the significance of observational 
constraints on the post-Einsteinian extensions. 



VI. PERIHELION PRECESSION 



We now study the influence of post-Einsteinian extensions of GR on the bound geodesies followed by planets. We 
discuss in particular the perihelion precession, an observable well known in tests of gravitation theory 1] . As for other 
observables studied in this paper, we focus our attention on anomalies induced by variations of the metric around the 
standard GR situation. 

As the non linearity of gravitation theory plays a key role in the problem of perihelion precession, we present quite 
detailed calculations, starting from the Hamilton- Jacobi equation written, for a static isotropic metric, in terms of 
time i, azimuth angle and a phase shift 5,- {E and J being the constants of motion H26|l ') 

S = -Et + Jip + Sr 

{drSrf - -grr ~ A ~ (39) 
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To analyse bound orbits, it is convenient to use Schwartzschild coordinates and rewrite the quantities (|39|l in terms 
of corresponding metric components 




Sr-TJl J-f^W^-5oo(^+^^l§ (40) 



^-.goo(^ + s2) Jc'J -goon' 

We recall that u is the inverse of the radial coordinate f and that the prime symbol used below denotes derivation 
with respect to u. 

One obtains in particular the orbital period, that is to say the duration between two successive identical values of 

(fi, as 

E f^'^ dip , , 

T = / ^ 41 

For circular orbits with a constant radial velocity, t and ip are proportional to each other, and the orbital period only 
depends on the metric component goo- Equations 140|l and (|41(l thus lead to the third Kepler law {u being the inverse 
radius of the circular orbit) 

T=— — , +u^^-^—^~—goo 42 

c Vffoo/ J J^c^ 900 900 

Hence, the orbital period may be written as the sum of a standard GR expression and of a post-Einsteinian anomaly 

T=[Tl,+5T , [Tl,=2.n-H-i , ^ = --^ (43) 

These relations entail that third Kepler law depends on variations of goo only. In other words, bounds resulting from 
third Kepler law constrain a particular linear combination of the two independent possible variations. We will come 
back on this discussion later. 

We now study perihelion precessions which can be determined from the variation A(/5 of the azimuth angle between 
two identical values of the radial coordinate f, the elapsed time being nearly but not exactly one orbital period. This 
quantity may thus be obtained from the phase shift Sr (|40|l evaluated over this elapsed time 



goo V J2c2 ^'"^^ J2 ' '^2 

djASr = 2t: + Aip (44) 

We now evaluate the variation of these quantities for a perturbation around GR, the constants of motion E, J being 
kept fixed. We obtain the following anomalies for the phase shift and precession angle, at first order in the metric 
perturbation. 



ASr = [AS'rlst + ^^^r , Aip = [Ap]^^ + SAp 

SASr ^ -\j [ffOolst ^9rrd [Sr]^, + j [g^f U^ '^ 

5Ap} ^ djSASr (45) 



For the nearly circular orbits followed by planets, the anomalous precession SAip may be expanded in the orbit 
eccentricity e. This expansion is presented in IdI and we give here its result for the post-Einsteinian anomaly of the 
perihelion precession 

^ ( X. + f!^ ( 5-g';^ - ^\ \ (46) 



[gool, 4 \ - 2[g, 



OOJst 
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Inserting explicit expressions for metric perturbations H19I20|I in terms of the potentials and 5^p (with S^n 
considered to be much smaller than 5$p), we also obtain 

^ ^ . {u5^p) - ^ + ^ [{u'S^'i^ + ^^^'i'^) - j (47) 

At the limit of nearly circular orbits e 0, the precession anomaly H47|l takes a simple form depending only on the 
second derivatives of the functions S<^n and uS<^p 

C-^) ^ {2.uS^P~S^,)"^ (48) 

We may recover the well known expression of the anomalous perihelion precession in the PPN formalism 1] by 
inserting in H48I) the particular expressions (|C3|I of S'^jy and (5$p for the PPN metric. Keeping only the leading term, 
we obtain the anomalous perihelion precession in terms of anomalous parts of Eddington parameters /3 and 7 



['^ ) ^ {2ij-l)-if3-l))Ku (49) 

The more general expression (|48f) may thus be thought of as stating that the functions (5$'^/2k^ and {uS^p)" /2k 
promote the anomalous PPN parameters /3 — 1 and 7 — 1 from the status of constants to that of space dependent 
functions. However this identification can only be a first rough description of the more general expressions obtained 
throughout the paper. In the next section, we use these more general expressions to discuss how planetary tests fix 
bounds on potential post-Einsteinian extensions of GR. 

VII. OBSERVATIONAL CONSTRAINTS 

We now use the expressions obtained in previous sections to draw constraints on post-Einsteinian extensions of 
GR from observations or experiments performed in the solar system. Any potential anomaly exhibits a particular 
dependence on the two functions of the radius coordinate used to parametrize the extensions. This freedom will 
be the key point allowing us to accomodate the Pioneer anomaly besides the other gravity tests performed in the 
solar system. The two independent functions can be chosen as the anomalous metric components {Sgoo^^grr), the 
potentials [6'^n,S'^p) or the Einstein curvatures [SEq, SEl/Y Relations between these representations have been 
given in sect. 3, so that they can be used equivalently for the calculations. For discussing the significance of the tests, 
we will emphasize the parametrizaton in terms of anomalous Einstein curvatures since they describe in the clearest 
manner the deviation from GR. 

We first discuss the Pioneer-like anomalous acceleration (|38|) calculated for probes with an eccentric motion in the 
outer solar system. At leading order, this expression can be linearized since non linear corrections would involve a 
small Newton potential at large heliocentric distances. We obtained in this manner a simplified expression 

Sa 1„, f dr\^ 1 ,^ (fr 

?oo ' ' 



^ drS^N~2dr{6^N-S^p)i — ] - (25$^ - <5l>p) — (50) 

Now, the first contribution, which depends only on S^n, can simply be interpreted as resulting from a long range 
modification of the Newton force law, and it cannot explain the Pioneer anomaly without being noticed on planetary 
tests 0,1^. Then, the last term proportional to the acceleration of the probes can be ignored, since the latter keep 
practically constant velocities over the observation range. Finally, the term (5$jv can be neglected in front of 5$p in 
the parenthesis appearing in the second term 16, 17]. We thus obtain 

Sa ~ 2vldrS^p , VP ~ 1.2 x IQ-^ms^i (51) 

where vp is the velocity of Pioneer probes. If the result is identified with the anomalous acceleration recorded on 
Pioneer probes, we get an information about the function 6^p at large distances 

6a~~8x 10"^°ms-2 , c^u^duS^p = -c^drS^p ~ 0.25ms~2 (52) 
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Pioneer data indicate that the anomalous acceleration remains practically constant over a large distance range 
20 — 70 AU of heliocentric distances in the outer part of the solar system. This suggests that the potential (5$p 
has a linear dependence in the radial coordinate f or, equivalently, varies as the inverse of u over the same distance 
range. This behaviour appears to favor a simple post-Einsteinian model Il7j characterized by a single additional 
parameter with respect to GR, the latter being fixed by ifS^ . 

(^*^).nodei--^ ' CpM-0.25ms-2 (53) 

It should be emphasized that a similar component introduced in the potential is constrained by planetary tests 
to remain vanishingly small 

('5^-),„odei--^ . < 10--ms- (54) 

This estimation given in should be refined after the discussions of the present paper but the difference in the 
orders of magnitude is large enough to ensure that 5^n is negligible in front of S^p in the outer solar system. We 
may then derive from H21|l the expressions of the anomalous Einstein curvatures 

fXT?o\ 'KpMu 2CpMu 

As already emphasized, this is the clearest manner to characterize the post-Einsteinian extension of GR studied in 
the present paper. In contrast to GR, the extended theory corresponds to Einstein curvatures with non null values far 
from the gravity sources. This feature is not a complete novelty since the PPN metric also shows Einstein curvatures 
differing from their standard null values in free space surrounding the gravity sources. But PPN expressions (|G2|I 
exhibit peculiarities which entail that they cannot account for the Pioneer anomaly. First, they lead to a non null 
value of SE"!^ at first order in k but to 6Eq vanishing at the same order: SEq is negligible with respect to 5E1^ in 
the PPN metric whereas the two quantities have comparable magnitudes in 1)55(1 . In other words, the PPN Ansatz 
explores a particular linear combination of the two possible extensions which is not the linear combination producing a 
Pioneer-like anomaly. Then, the anomalous Einstein curvatures obtained in the PPN metric scale as for the leading 
term, which is again clearly distinct from those associated with the model (|53() which scale as u. In other words, the 
PPN metric which is defined as an expansion of the metric in terms of Newton potential can only produce short-range 
effects whereas the Pioneer anomaly is rather a long-range modification of GR. Note that such a long-range effect can 
only be obtained with an infrared modification of the running couplings [l6l Il7l | . 

The preceding discussions refer to the values of the two potentials 6^]^ and 6^p at large heliocentric distances, 
typically in the range 20 — 70 AU where the Pioneer anomaly has been recorded. But the presence of the same 
potentials might also be detected in the inner part of the solar system where very accurate measurements have been 
performed. This question was discussed in the end of ^3 for the case of deflection experiments which can be analyzed 
within a linearized framework. Here, we consider planetary tests bearing on third Kepler law and perihelion precession 
Ub]. As already discussed, tests of the third Kepler law H43|) depend on variations of the temporal metric component 
§00 only. At leading order, these variations amount to a modification of the Newton potential S^n which must remain 
very small in order to reproduce the agreement of GR with planetary measurements [Tol Isif . This leads to the bound 
()54|l written above for the function (5$jv. 

The discussion for the bounds drawn from perihelion precessions involves not only the modification of the Newton 
potential 5$jv but also the second potential (5$p. In fact the constraints on (5l>jv prevent it to significantly contribute 
to the perihelion precession anomaly, so that H47|l can be translated into 



27r 



2 

{u5^p)" -f ^ (^u{u5^p)"')' + u {u5^p)" - 2u5^'p) (56) 



We have kept the leading term which subsists at the limit of a circular orbit as well as the sub-leading terms which 
depend on the squared eccentricity. Noticing that the dominant term in H56|l vanishes for the simple model l|53|) 
suggested by the Pioneer anomaly, we deduce that the sub-leading one thus determines the perihelion precession 
anomaly. To give an idea of the orders of magnitude, we write H56|l for the simple model H53() . considering the Mars 
orbit (radius tm 1.5AU, eccentricity e^j ~ -09) and assuming that the parameter Cp has the value (|53|l needed to 
fit the Pioneer anomaly 

^^-^^^""^^^1.2x10-^ (57) 
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This value is of comparable order of magnitude than the accuracy on Mars perihelion precession (see the bound 
deduced in from the data of jl^)- This means that perihelion precession anomalies might be used to detect the 
effect of the second potential in a dedicated analysis starting from the new model . 

This analysis would have to look for values of the new parameters added with respect to a PPN model: the first of 
these parameters is but the variation of Cp versus distance should also be taken into account. To make this point 
clear, let us consider a model differing from H53|l through the addition of other power law terms 

('5*^).„odcl = -22 



The term n = in this sum reproduces 1)5311 and leads to a perihelion precession anomaly proportional to the squared 
eccentricity. The other terms n ^ lead to relatively larger contributions to the perihelion precession anomaly 

[^J .o.e> ^ - ? 0^ I" + ^ + 1) + 2) j (59) 

To fix ideas, let us calculate the effect for the Mars orbit, taking into account the next term n = 1 besides the term 
n = 0, 

\ / model \ / L,p 

As already alluded to, the perihelion precession anomaly may thus be appreciably affected (for values of e of the 

order of 2|i) though the potential remains close to H53|) (because e <C 1). In other words, the perihelion precession 
anomaly could be used as a sensitive probe of the parameter Cp and of its relative variation e. We stress once again 
that these values are evaluated at the radius of Mars orbit, if the perihelion precession of Mars is studied, and not at 
the heliocentric distances of Pioneer probes. Dedicated planetary tests could therefore bring to our knowledge new 
informations on the potential deviations from GR, the aim being to cover the intermediate distance range between 
the close vicinity of the Sun, probed in deflection experiments 17], and the outer solar system, probed by Pioneer 
spacecrafts. 



VIII. CONCLUSION 



In this paper, we have presented a new kind of extensions of GR which allows one to explore potential variations of 
gravitational theory. These extensions, which are motivated by the radiative corrections arising from the coupling of 
gravity field to the other fields, have been explored from a phenomenological point of view. Their ability to account 
for the Pioneer anomaly while remaining compatible with other gravity tests performed in the solar system, which 
had been proved within a linearized treatment of gravitation fields 16, 17], has been confirmed here without this 
restriction. 

The new framework generalizes the PPN Ansatz jj as well as the idea of a long-range modification of the gravity 
force law The new framework is characterized by two potentials which are functions of the radius: the first one can 
be identified with a modification of the gravity force law while the second one replaces the PPN constant 7. In fact, 
the first potential (5$ n can only slightly deviate from its standard form, as a consequence of the good agreement of 
planetary tests with GR. The key ingredient allowing one to account for the Pioneer anomaly is the second potential 
(5<lp with a non standard long range dependence, approximately a linear variation of the second potential versus the 
radius. We have shown that the PPN metric tensor was a particular extension of GR which did not show the ability 
to account for the Pioneer anomaly, in contrast to the general post-Einsteinian metric. 

The second potential S<tp, needed to reproduce a Pioneer- like anomaly, also leads to other observable anomalies, 
in particular in time delay, deflection or Do pplc r tracking experiments on electromagnetic sources occulted by the 
Sun. Such anomalies have been analyzed in |l7| and the results of this analysis are unchanged because a linearized 
treatment is sufhcient to this aim. We may recall here that crucial indications could be obtained through a reanalysis 
of existing data, in particular those of the Cassini experiment 139]. On a longer term, they can be drawn from high 
accuracy Eddington tests (see for example the LATOR project [53) or global mapping of deflection over the sky (see 
for example the GAIA project [13] )■ 
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In the present paper, we have studied the effect of the second potential S^p on the anomalous perihelion precession 
of planets, in particular Mars. We have shown that dedicated analysis of planetary data could also bring new 
information of interest on the variation of S^p at planetary ranges. Another possible test of the new framework could 
be to confront its detailed predictions with the sets of Doppler tracking data recorded during the flight of Pioneer 
10/11 probes towards the frontiers of the solar system. Such a study might provide interesting clues, in particular 
through an inspection of the variation of the anomaly during and after the flyby sequences which brought the probes 
to their escape trajectories. Another attractive idea would be to confront the predictions of the new framework with 
the results of a specifically designed mission . 



APPENDIX A: VICINITY OF EINSTEIN THEORY 

In this appendix, we recall the general definitions of curvatures, with the conventions of [4^ . and we write Bianchi 
identities obeyed by these curvatures. We then deduce simplified relations valid for gravity theories lying in the 
vicinity of GR. 

We first reproduce the definitions of Christoffel symbols F^^, Riemami (Rxfu^p), Ricci [R^u) and scalar {R) curva- 
tures 

A 9^^ 
Tp,. = —{d^^gup + dyg^^p-dpg^u) 

— c\ pA o pA I pA po" pA per 

^pyp — t^i/J- pp ^p-L pi/ I -L pp ^ pa pu 

Rpv = Rp\„, R = Rp (Al) 

We also recall the definition of Einstein curvatures Ep^, 

1 

r 

and Weyl curvatures Wxp^p 



Ep, = Rp, --R, E = -R (A2) 



Wxpup = Rxpup — - {gxvRpp + gppRxu — gxpRpv — gpvRxp) 

+^ig\i^gt^p- gxpgtj^iy) R (A3) 

Riemann curvatures express the commutators of translations of vector fields in a curved space j5(|. They satisfy 
Bianchi identities, that is also Jacobi identities for the corresponding algebra, 

^crRpp,^ + ^yRpap + ^pRp^a = 

^xRpp. = Vpi?^, - VyRpp , VxR^p = l^pR (A4) 

These identities entail in particular that Einstein curvature tensor is transverse (V^ is the covariant derivative asso- 
ciated with the metric) 

S/xE^ EE dxEl + r^^i?^ - VP^^E^ = (A5) 
They also lead to relations between covariant derivatives of the Weyl curvatures and Einstein tensor 

VaW^p'.p = ^V, (epp - \gppE\ - iv. Up, - ^gp^E 



V-VaT^X - l^'Epp + ^{\/p\/p-gpp\/')E-^W;^^,E: 

+E;E^p ~ \eEpp - gpp - ^) (A6) 

The latter identities have a form well suited for an expansion of gravity theory around GR. In the latter theory 
indeed, Einstein curvatures vanish everywhere except on gravity sources, that is the Sun in the solar system. This 
entails that identity (jA5l) thus has a simplified form 

[Va],, SE^ e. dx5El + [Vi;\ 5EP^ - [r^ J 6E^^ = (A7) 
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The simplification is even more important for the identity (|A6|) with the disappearance of quadratic forms of Einstein 
curvatures 



(A8) 



Note that these equations are not vahd on the gravitation sources. They can be simphfied one step further by using 
the fact that potentials or curvatures have a relatively small amplitude in the solar system fsee lB)) . 

APPENDIX B: APPROXIMATIONS IN THE SOLAR SYSTEM 

In this appendix, we write down approximated expressions by using two assumptions which are valid in the solar 
system: the extended gravitation theory lies in the vicinity of GR and curvatures have a relatively small amplitude. 

The first assumption allows one to use expressions obtained in^ The second assumption then leads to a substitution 
of covariant derivatives by ordinary derivatives. The transversality of Einstein curvatures (jA7|) is thus simply read as 
the ordinary Minkowskian expression 

dx5El = (Bl) 

Using also the other Bianchi identities l|A8|l . the transverse Einstein tensor may then be decomposed using projectors 
(^0) and (ttI) 

5E^, = <5(i;°)^,+^(£;i)p, (B2) 
5{E'')^, = {^'')^,,xp6E^P , 5{E^)p, = {^^)p,xp5E^'' 

[T^ jpvXp — ^ ^ 

\" )pvXp — 2 J i^nu — 'Itiv Q2 

Equations IjASp relating Einstein curvatures to Weyl and scalar curvatutres then take a simple form 

5Ep^{x) = -{■k°)p_i,qqA5w{x) - {tt^)p.i,ooSR{x) 

A = d^ + -dr (B3) 
r 

This reproduces the decomposition on the sectors of traceless and traced tensors which was used in [l6l IT^ . We recall 
that the indices refer to the conformal weight of the components, for traceless tensors and 1 for traced ones [35l| . 

The modification of gravitation equations may also be characterized in the momentum domain by introducing scale 
dependent running coupling gravitation constants and which differ in the two sectors [l^ [T^l ■ The perturbed 
Einstein gravitation equations H12|) may indeed be rewritten (compare with sect. 2 in [IT)) 

/(-) - / ' - - ^ 

rr. r,i . 0^ 8TT6G°[k]M . ^. 8TT6G^[k]M 

SEp^[k] = (7rO)^^oo ^j2^— + ('^ )m'>oo j^^" (^4) 

The information contained in the anomalies SG'^ and 6G^ of the running gravitation constants is completely equivalent 
to the spatial dependence of the curvature components, for instance Sw and 6R, 

ASwix) ~ SR{x)~ (B5) 

c 

The same information may equivalently be presented in terms of the modifications of Einstein components 

_o 2 1 2S& + SG' 8nM 
6E„ ~ --A(5w - ~ — 



2^11,^^ 2 ^ 1 6G^ - 6G° 8nM 

5E; ~ -dr--{A5w - 5R) -dr-r — B6 

r 3 A r A 3 
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or, alternatively, in terms of two potentials 5^n and 5^p (see eas. (|18ll : Schwartzschild coordinates are used) 

AJ$jv - 5 T- ' A = d^ + -dr 

id' r 

-r,- 2{SG" ~ SG^) 47tM 

AS^P 0, ^ (B7) 

In particular, potentials varying linearly with the radius [l^ IT^ or, equivalently, Einstein curvatures varying as 
the inverse of the radius correspond to anomalous couplings 5G^[k] and (5G^[A:] scaling as 1/k^ where k is the spatial 
wavevector. 

APPENDIX C: THE PPN ANSATZ 

In this appendix, we write down the PPN ansatz as a particular case of the more general framework presented in 
the rest of the paper. 

The PPN metric is written in terms of isotropic coordinates with 



300 = 1 + 20 + 2/30^ + . . . , -Ku 

grr = -l+27(/.- ^(502 + ... (CI) 

The Eddington PPN parameters (3, 7 and 5 are equal to 1 in general relativity and close to 1 in its vicinity. The 
symbol . . . has been introduced to recall that these expressions are but the first terms in an expansion in terms of the 
Newton gravitation constant Gn or, equivalently, of the Newton potential 0. 

It is instructive to write down the expansion of the anomalous Einstein curvatures 

5E;~ 2(7 - 1)ku3 + (4(^_1)_ 12(7-1) + 3(5 -l))K2ii'* + ... (C2) 

This expansion shows that 5E'^ vanishes at leading order in Gjv, whereas 5E^ does not. In other words, the PPN 
metric is a particular extension of GR which affects (5i?^ more than 5E^. 

It is also worth writing the anomalous parts of the two potentials 5$ n and (5$ p for the PPN metric 

5#(^^^) ~ (/3 - I).c2s2 (C3) 

<5#r"^ - 5#r"^ -(7 - i).u + g(7 1)) 

This shows that (5$jv and S'^p vary quadratically with the inverse radius. In other words, the PPN metric is a 
particular extension of GR with anomalies present at short ranges (in the inner solar system) rather than at long 
ranges (in the outer solar system). 

APPENDIX D: PRECESSION ANOMALY 

In GR, the gravitation potential takes its standard Newtonian form. We may then solve the Kepler problem by 
introducing a polynomial \V]^^ of third degree in the inverse radial coordinate (Schwartzschild coordinates are used 
in this appendix), with coefhcients determined by the two invariants J and E 

E^ , fm^c^ 



[^]st (S) ^ ^ - [.9oo].t ( ^ + ^ ) (Dl) 
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We may also write it in terms of the Keplerian elements, namely the orbit mean inverse radius p and the orbit 
eccentricity e, 

2a 



[VI, [u) = [ eV -{u-pY + —{u~ pY 



J2 

Hamilton- Jacobi equations ()39|l lead to the following solution 



A ^ 3 , \ 1 , 
P ( 1 + 2 iaooU.P j =^2 ~P' 

nHe^-l~-2a)p^^^-^ (D2) 



which allows one to recover the perihelion precession for Einstein geodesies 



2^ + [A^],, = <f d i^l, =±j (D4) 



1 f^^ (1 — 2eo'sin(y5) ^cosipdip 



cos2(/7 + j5^(l — 6(7 — -^1 — 2ecrsin(p)2 



We come now to the post-Einsteinian case, and expand the generating function for precessions, namely AS'j. H44|l . 
in the vicinity of GR metric (|45|l . This function depends on the constants of motion and metric components through 
two generic integrals 

5^Sr ^ ^^IsilMsJ-grA-^^I^i^^} (D5) 

I IJc [gQQ\^, 



Is{f} ^ j f{u)'^^ , lM}^ffiu)dH 



Corrections brought by the Newtonian potential are of order | [goo]ltPl ^ may be calculated in a perturbative 
expansion. Integrals in equations (|D5|I are thus expanded around expressions for Keplerian orbits {il — 1,(t = 0). 
Noting that / J^c^p^ is of order —I (see (ID2II '). one rewrites expressions in (jD5|) in terms of p and e 



P' [- ,- 1/ +e- -4) , a~-^[goottP 
' p{\^ esin(/3 + — e^sin^(p^ , d \^\, ~ (I + (t(3 + esiniys)) dip 



Ik: 



/{/}=/ /(p(l + esin^))d^ , /s{/}~/{p-p)2_eV]4} 
Jo 

74/} (1 + 3a)/{/} + ^/{((« - pfm (D6) 

In fact, confrontation with observations shows that the two potential perturbations do not occur at the same level, 
but that corrections to the first potential occur at one order higher than to the second potential. Taking this remark 
into account, and using relations (|D2|I . we obtain from ljD5(l and (|D6|1 the precession anomaly at leading order 

,5A^ ^ -\{\ + ja;)/s{%.} + -^(1 - 

Jdj ~ -2 (p9p + (1 - e2)5e2) (D7) 
For a nearly circular orbit, expressions (jD6p may be expanded in the orbit eccentricity e 



(1 - J9j)/{/} ^ (^(uV)"b) + 



2^2 

(D8) 
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One then obtains the perihelion precession anomaly at leading order 



SA(p 



-TT Sgrr - 




{ 



2[5oo]: 



St 



) 



) 



(D9) 



Functions of the radial coordinate have been evaluated at 1/p. 



REFERENCES 



[1] Will C. M., Theory and experiment in gravitational pliy.sics (1993) Cambridge U. P.; Living Rev. Rel. 4 (2001) 4. 
[2] Fischbach E. and Talmadge C, The Search for Non Newtonian Gravity (1998) Springer Verlag. 

[3] Perlmutter S. et al , Astrophys. J. 517 (1999) 565; Perlmutter S., Turner M. S. and White M., Phys. Rev. Lett. 83 (1999) 
670. 

[4] Aguirre A. et al , Class. Quantum Grav. 18 (2001) R223. 

[5] Ricss A. G. et al , Astron. J. 116 (1998) 1009. 

[6] Sanders R. H. and McGaugh S. S., Annu. Rev. Astron. Astrophys. 40 (2002) 263. 

[7] Luc A., Scoccimari R. and Starkman G., Phys. Rev. D 69 (2004) 044005. 

[8] Carroll S.M. et al , Phys. Rev. D 70 (2004) 043528. 

[9] Anderson J.D. et al , Phys. Rev. Lett. 81 (1998) 2858. 
[10] Anderson J.D. et al , Phys. Rev. D 65 (2002) 082004. 
[11] Anderson J.D. et al , Mod. Phys. Lett. A17 (2003) 875. 
[12] Nioto M.M. and Turyshov S.G., Class. Quantum Grav. 21 (2004) 4005. 

[13] Turyshcv S.G., Nieto M.M. and Anderson J.D., 35th COSPAR Scientific Assembly (2004) Preprint gr-qc/0409117. 

[14] Bertolami O. and Paramos J., Class. Quantum Grav. 21 (2004) 3309; see also Preprint astro-ph/0408216 and Preprint 

gr-qc/0411020. 
[15] Dittus H. et al , Preprint gr-qc/0506139. 

[16] Jaekel M.-T. and Reynaud S., Mod. Phys. Lett. A20 (2005) 1047. 

[17] Jaekel M.-T. and Reynaud S., Class. Quant. Grav. 22 (2005) 2135. 

[18] Utiyama R. and Do Witt B., J, Math. Phys. 3 (1962) 608. 

[19] Deser S. and van Nieuwenhuizen P., Phys. Rev. D 10 (1974) 401. 

[20] Capper D.M., DufT M.J. and Halpern L., Phys. Rev. D 10 (1974) 461. 

[21] t'Hooft G. and Veltman M., Ann. Inst. H. Poincare A 20 (1974) 69. 

[22] Stelle K.S., Phys. Rev. D 16 (1977) 953; Gen. Rel. Grav. 9 (1978) 353. 

[23] Fradkin E. S. and Tseytlin A. A., Nucl. Phys. B 201 (1982) 469. 

[24] Lauschcr O. and Router M., Class. Quantum Grav. 19 (2002) 483. 

[25] Reuter M. and Weyer H., Phys. Rev. D 70 (2004) 124028. 

[26] Bonanno A. and Reuter M., J. High. En. Phys. 0502 (2005) 035. 

[27] Adelberger E. G., Heckel B. R. and Nelson A. E., Ann. Rev. Nucl. PaH. Sci. 53 (2003) 77, also in Preprin£ hep-ph/0307284. 

[28] WilUams J. G., Newhall X. X. and Dickey J. O., Phys. Rev. D53 (1996) 6730. 

[29] Hellings R. W. et al, Phys. Rev. Lett. 51 (1983) 1609. 

[30] Anderson J. D. et al , Astrophys. J. 459 (1996) 365. 

[31] Jaekel M.-T. and Reynaud S., Int. J. Mod. Phys. A20 (2005) 2294. 

[32] Damour T., Class. Quant. Grav. 13 (1996) A33; Damour T., Piazza F. and Veneziano G., Phys. Rev. D 66 (2002) 046007. 
[33] Overduin J. M., Phys. Rev. D 62 (2000) 102001. 

[34] Hilbert D., Nachrichten von der Gesellshaft der Wissenshaften zu Gottingen (1915) 395; Einstein A., Sitzungsherichte der 

Preussischen Akademie der Wissenschaften zu Berlin (1915) 844. 
[35] Jaekel M. T. and Reynaud S., Annalen der Physik 4 (1995) 68. 
[36] Talmadge C. et al . Phys. Rev. Lett. 61 (1988) 1159. 
[37] Shapiro I. I., Rev. Mod. Phys. 71 (1999) S41. 
[38] less L. et al , Class. Quant. Grav. 16 (1999) 1487. 
[39] Bertotti B., less L. and Tortora P., Nature 425 (2003) 374. 
[40] Shapiro S. S. et al , Phys. Rev. Lett. 92 (2004) 121101. 

[41] Shapiro I. I., in general relativity and Gravitation 1989 N. Ashby, D. F. Bartlett and W. Wyss eds (Cambridge UP, 1990) 
313. 

[42] Nordtvcdt K., Phys. Rev. 170 (1968) 1186. 

[43] Landau L. D. and Lifschitz E. M., The classical theory of fields (1975) Butterworth. 
[44] Misner C. W., Thorne K. S. and Wheeler J. A., Gravitation (1972) Freeman. 



19 



[45] Eddington A. S., The mathematical theory of relativity (1957) Cambridge U. P.. 

[46] Robertson H. P. in Space age astronomy (1962) Academic Press. 

[47] Ross D. H. and Schiff L. I., Phys. Rev. 141 (1966) 1215. 

[48] Nordtvedt K., Phys. Rev. 169 (1968) 1014 & 1017. 

[49] Will C. M. and Nordtvedt K., Astrophys. J. 177 (1972) 757. 

[50] Nordtvedt K. and WiU C. M., Astrophys. J. 177 (1972) 775. 

[51] Reference 13, chap. 23. 

[52] Coy J., Fischbach E., Hellings R.W., Talmadge C. and Standish E.M., private communication (2003). 
[53] Turyshev S.G. et al , in Trends in Space Science and Cosmic Vision 2020 to be published (ESA, 2005), also in Preprint 
gr-qc/0506104. 

[54] Vecchiato A., Lattanzi M.G. and Bucciarelli B., Astron. Astrophys. 399 (2003) 337. 

[55] Dittus H. et al , in Trends in Space Science and Cosmic Vision 2020 to be published (ESA, 2005), also in Preprint 

gr-qc/0506139. 
[56] Reference H, chap. 15. 



